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Abstract. The finite lattice method of series expansion is generahsed to the qf-state 
Potts model on the simple cubic lattice. It is found that the computational effort grows 
exponentially with the square of the number of series terms obtained, unlike two- 
dimensional lattices where the computational requirements grow exponentially with 
the number of terms. For the Ising [q = 2) case we have extended low-temperature 
series for the partition functions, magnetisation and zero-field susceptibility to u^^ 
from The high-temperature series for the zero-field partition function is extended 
from v^^ to v'^'^. Subsequent analysis gives critical exponents in agreement with those 
from field theory. 
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1. Introduction 



This is the first of a series of papers describing the application of the finite lattice 
method of series expansion to various cases of the g-state Potts model. A major 
objective in this work is the study of the q = 3 case in three dimensions because of 
its relevance to quantum chromodynamics. However in the course of this project we 
have developed extremely powerful techniques of series expansion. This has made it 
possible to obtain considerable extensions to many Potts model series and so we have 
been able to investigate a large range of lattice statistics problems. 

In this paper we describe the application of the finite lattice method of series 
expansions to the derivation of high- and low-temperature expansions for the free 
energy of the Potts model on the simple cubic lattice. We present series for the 
Ising model which is the q — 2 case of the Potts model. We have been able to 
extend the low-temperature series for the zero-field partition function, magnetisation 
and susceptibilty from order tt^° (Sykes et al 1973) or tt^^ (M.F. Sykes, unpublished) 
to u^^. The zero- field high-temperature series for the partition function has been 
extended from v^^ (Sykes et al 1972) to v^^. The variables u = exp(-4J/A;T) and 
V = tanh( J/fcT) are the usual low- and high-temperature model expansion variables for 
the Ising model. (A recent paper by Bhanot et al (1992) calculates low-temperature 
series for the simple cubic Ising model using a variant of the finite lattice method that 
is apparently less efficient than ours. They obtain series for the magnetisation to u'^^ 
and the internal energy to order w^^, thus duplicating the work of Sykes et al (1973) 
and Sykes (unpublished) respectively.) 

The extension of these series is particularly pertinent as their behaviour is quite 
difficult to determine. While the high-temperature series are well-behaved, with 
uncertainty as to critical exponents being confined to the fourth, or at the worst 
third, decimal place, the situation at low temperatures is far less satisfactory. For 
many years there was considerable controversy as to the symmetry or otherwise of 
the critical exponents above and below the critical temperature. That this is now 
less controversial is not based on good numerical evidence however, but rather on a 
strengthening belief in universality based on a better understanding of the connection 
between field theory and critical phenomena. Nickel (1991) has also recently argued 
for the extension of perturbation series, in order to remove any lingering doubts about 
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the equivalence of the 0^ field theory, and the lattice models. 

The new low-temperature series we have obtained here permit an improved 
analysis, though the quality of the numerical estimates is still inferior to that of the 
high-tcmpcrature series. Nevertheless, by biasing our analysis with the value of the 
critical temperature, obtained from the high-temperature series, and by assuming a 
confluent exponent around | (the exact value is unimportant), exponent estimates 
are obtained that are consistent with the best field theory estimates. We have also 
extended the high-temperature zero-field free energy series, which allows us to estimate 
a considerably more accurately than has been possible previously. 

The finite lattice method of series expansions for lattice statistics problems has 
been applied to a number of different systems on a range of two-dimensional lattices. 
The initial application was in calculating the high-temperature expansion of the three- 
state Potts model (de Neef, 1975; de Neef and Enting, 1977). The expansion of 
the limit of chromatic polynomials (Kim and Enting, 1979) was formally a high- 
temperature expansion and the enumeration of self-avoiding polygons (Enting, 1980a; 
Enting and Guttmann, 1985; Guttmann and Enting, 1988a) is closely related to high- 
temperature expansions. However most of the subsequent applications of the method 
(based on the work of Enting, 1978a) have been in the derivation of low-temperature 
or high-field expansions, generally on the square lattice (e.g. Enting, 1980b; Adler et 
al 1983). Series for triangular lattice systems have been obtained by regarding them 
as square lattices with additional interactions (Enting, 1980c; Enting and Wu, 1982). 
Specific combinatorial expressions for triangular lattice systems are known (Enting, 
1980d, 1987a). However these involve the additional complication of calculating 
the partition functions of hexagonal finite lattices. The only application using full 
triangular symmetry has been in the enumeration of polygons on the triangular 
lattice (Enting and Guttmann, 1992). The finite lattice method has also been applied 
to Potts models on the checkerboard lattice (Enting, 1987b) and polygons on the 
honeycomb lattice (Enting and Guttmann, 1989). These last two cases were treated 
as modifications of the square lattice. 

Enting (1978b) quoted some of the combinatorial results required for the finite 
lattice method on the simple cubic lattice but until now these results have not been 
applied. In this paper we describe the formalism for obtaining either high-temperature 
or low-temperature expansions for the Potts model on the simple cubic lattice. We 

4 



present and analyse specific results for the Ising model case. 

The outline of the remainder of the paper is as follows. Section 2 analyses the 
combinatorial aspects of applying the finite lattice method on the simple cubic lattice. 

Section 3 describes the expansions that we have calculated. Section 4 describes an 
analysis of these series. 

2. Finite Lattice Expansions on the Simple Cubic Lattice 

We begin by defining the notation that we shall use throughout this series of 
papers. The Potts model is defined on a lattice in terms of a 'spin' variable, Sj at 
each site, j, taking integer values from to q — I. There is an energy difference, 
AE, between aligned and non-aligned states that interact. Generally we consider 
interactions confined to pairs of nearest-neighbour sites, i.e. those joined by the bonds 
of the lattice. Sites aligned in the '0' direction differ by a field energy H from sites with 
other alignments. Thus with S{m,n) — 1 of m — n and otherwise, the Hamiltonian 
is written 



where the first sum is over all pairs of interacting sites and the second sum is over all 
sites. 

The low-temperature expansion variable is 2; = exp(— /3A£^) where /5 = 1/kT. For 
the field dependence, we use the expansion variables /i = exp(— and x = 1— /i. The 
general Potts model high-temperature expansion variable is = {1 — z)/{l + {q — l)z). 

For the Ising model (i.e. g = 2), only even powers of z occur and the natural 
low-temperature expansion variable is tt = 2^ = exp(— 4/3J) and the high-temperature 
expansion variable can be written as v = tanh(/3J) with 2J = AE. For zero-field 
on bi-partite lattices, only even powers of v have non-zero coeffcients in the high- 
temperature Ising model expansion. 

The basic formulation of the finite lattice method approximates the partition 
function per site, Z, as 



(1) 



(id) 



Z = lim Z-, 
irHoo 



i/|r| 



n Z^^""^ 



(2a) 
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where F denotes a lattice that becomes arbitrarily large in all directions and |r| denotes 
the number of sites in F. Here ^ is a set of finite lattices, a, with A closed under the 
operation of intersection of finite lattices. For the simple cubic lattice, this general 
relation has the specific form: 

Z = hm Z^/^; ^ n ^^J"'''''^ (2^) 

AT— »oo 

[p,q,r]&A 

where Zpgr is the partition function of a cuboid of dimensions p x q x r sites. For 
low-temperature expansions, the Zpqr are to be evaluated with a surrounding layer 
of fully-ordered sites. For high-temperature expansions, the Zpq^ are to be evaluated 
with free boundary conditions. (It is also convenient to remove common factors as 
described below.) The weights W{p, q, r) depend on the set. A, over which the product 
is taken. In approximations (2a, 2b) an appropriate choice of weights will give Z as 
a series correct up to, but not including, the order of the first connected graph that 
will not fit into any of the cuboids of set A (Enting, 1978a). 

The method of Bhanot et al (1992) uses generalised helical boundary conditions 
imposed in a sequence of configurations so that they can ultimately remove the effect 
of unwanted graphs that occur under helical or periodic boundary conditions. 

Inspection of the low-temperature expansion of the Potts model shows that the 
limiting graphs are trees that do not double back in any direction: all planes drawn 
perpendicular to bonds of the lattice intersect such trees at most once. Such a tree can 
span a cuboid of size pxqxr with, p + q + r — 2 sites and p + q + r — 3 bonds in the tree 
and will give powers of 4(p + g + r) — 6 or more in the Potts model low-temperature 
variable, z. 

Therefore it is appropriate to take A as the set of cuboids whose spans, cr, satisfy 

a— p + q + r<s 

For this choice of A, the first incorrect low-temperature term is of order 4(s + l) — 6, i.e. 
the series is correct to order 4s — 3. We use the notation A{s) — {\p,q,r] : p+q+r < s} 
to denote the set of cuboids used. 

Inspection of the high-temperature expansion shows that, for general fields, the 
same tree graphs are limiting. However, for zero field, the limiting graphs are 
maximally-extended polygons. These are polygons that have at most two intersections 
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with any plane perpendicular to the lattice bonds. These are the three-dimensional 
generalisation of the convex polygons (Guttmann and Enting, 1988b) that were used 
to obtain corrections in the polygon enumeration by Guttmann and Enting (1988a). 
For cuboids of span s + 1 the maximally-extended polygons will have 2(s -|- 1) — 6 steps 
so that cuboids of span < s will give the zero-field high-temperature series correct to 
2s — 5. The combinatorial factors from Enting (1978b) give: 



W{d, w,i)^ r]{p - d) r]{q - w) r]{r - £) for [d, w, L] e A{s) (3) 

\p,q,r]eA{s) 



where 



This implies 



r/(0) = 1 (4a) 

77(1) = -2 (46) 

rj{2) = 1 (4c) 

r]{k) = otherwise. (Ad) 



W{d, «;,£) = 1 for d + w + e = s (5a) 

= -5 foTd + w+£^s-l (56) 

= 10 ioid + w + e = s-2 (5c) 

= -10 foid + w + e = s-3 {5d) 

= 5 iord + w + £^s-4: (5e) 

= -1 ioTd + w + e = s-5 (5/) 

= otherwise. (4^) 

In actual computation it is convenient to exploit the cubic symmetry and consider 
only d <w < £. We define B{s) = {[p, q,r] : p + q + r < s,p < q < r}. The expansion 
becomes 

n ^P.r'"'''^ (6a) 

\p,q,r]eB{s) 

with 

V{p,q,r)= W{7r(p,q,r)) iorp<q<r (66) 

■K{p,q,r) 



where the sum is over all distinct permutations, 7r{p, q, r) of the indices. 

The partition functions are constructed by using a transfer-matrix formalism to 
build up £ layers of size d x w. As in all of the most recent applications of the finite 
lattice method, we used the approach of building up the finite lattices one site at a 
time. The computational complexity of the calculation is determined by the value of 
d X w that is required and this is determined by the limiting span s. Exploiting the 
cubic symmetry allows the series to be generated using cuboids whose size is subject 
to: 

1 < cZ < (imax and d < w < [{s — d)/2\ and w<£<s — d — w 

If s = 3m + 1 then d^ax = and the maximum size of (i x w is m x m. If 
s = 3m — 1 then dmax = m — 1 and the largest size of (i x w is (m — 1) x m. (The 
case of s = 3m has d x w = m x m but this is of little interest because the transfer 
matrices will be the same size as is required for the case s = 3m + 1 which gives 4 
extra low-temperature series terms and 2 extra high-temperature terms.) 

The limiting factor in the finite lattice series computations is the size of the vector 
required to store the site configurations at the end of a partly-built lattice. In a 
departure from the earlier high-temperature calculations noted in the introduction, it 
has been desirable to use a site representation rather than a bond representation 
when calculating the Zpqr required for high-temperature expansions. Thus size 
considerations are the same for both high- and low-temperature expansions. For 
the expansions of g-state models on the simple cubic lattice the number of site 
configurations required is q'' with k = d x w. For the g-state Potts model (g > 3) 
the symmetry between the states can be used to remove redundant elements. There 
will be approximately g'^/(g — 1)! distinct vector elements. The precise number can be 
calculated using a recurrence relation. Since the present paper only considers g = 2 for 
which exactly 2^^ states are required, the general recurrence relation will be described 
in the second paper. For zero-field high-temperature expansions, there is no need to 
distinguish the '0' state and so for a cross-section of k sites, approximately q^ /q\ vector 
elements are needed. 

2 

Thus vectors of approximately q"^ /(g — 1)! elements give Potts series correct to 
^i2m+i -^y^iIq vectors of q'^i'^-'^) /(g — 1)! elements give Potts series to z^'^'^~'^ . Vectors 
of size 2"^^ elements give Ising series correct to tt^"* while vectors of size 2™*^™~^) 
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elements give Ising series correct to u^'^~^. For zero-field high-temperature expansions, 
vectors of approximately q'^^ /q\ elements give Potts series correct to v^'^~^ while 
vectors of q'^^'^~^^ /q\ elements give Potts series to v^'^~'^ . Vectors of size 2"^^"-^ 
elements give Ising series correct to while vectors of size 2"^ -m-i elements 

give Ising series correct to 

In contrast, for the square lattice, the cross section, /c, corresponds to the width 
of the lattice (expressed as a number of sites) and so vectors of size f» q^ /{q — 1)!, 
will give Potts model low-temperature series correct to 4/c + 3. Duality gives the high- 
temperature series to the same order. Table 1 lists the low-order values of m together 
with the number of series terms obtainable for models on the simple cubic lattice and 
the size of vector required for the Ising model. 

While Bhanot et al (1992) compared the combinatorial complexity the 
combinatorial complexity of their technique to the finite lattice method with periodic 
boundary conditions, they made no comparisons with our formalism (see Enting, 
1987a; Guttmann and Enting, 1990b) using fixed boundary conditions. Their note 
gives too few details of the characterstics of their cancellation procedure for us to 
make a general comparison. However from the specific example, it appears that they 
need 24 spins (i.e. 2^^ configurations) to reach tt^^, while, as shown in Table 1, we 
require only 16 spins. 

In principle, it is possible to consider extending the series by determining correction 
terms to the finite lattice expansion. Previous studies (Enting and Wu, 1982; 
Guttmann and Enting, 1988b) have obtained explicit expressions for the leading- 
order correction terms for various models. For the minimal spanning trees direct 
enumeration seems difficult. However examination of the series shows that M such 
trees will give a correction of 

(g - l)Mi,^-'[z'^-' + i(g - 2)(. - 2)z^'-'] 

to the low-temperature series for the g-state Potts model. For g > 3 the number 
M can be determined by noting the correction required at the corresponding s value 
for the less complex q — 2 (Ising model) case. For the high-temperature series, the 
zero-field correction from maximally extended-polygons is of the form K(q — l)v'^'^~^. 
Again the Ising case can be used to determine the correction for higher g. Since the 
smaller the g the larger the lattice that can be used, our highest-order Ising cases can 
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be used to provide the correction terms for future q = 3 studies, as used in paper III 
of this series (Enting and Guttmann 1993, in preparation). 



3. Expansions 

The finite lattice method obtains the high- and fow-temperature expansions for 
the partition function, Z. For actual calculations, specific normalisations must be 
chosen. For low-temperature expansions, the appropriate choice is to define the 
partition function such that the fully aligned (all sites '0') state is taken as having 
zero energy, since the low-temperature series is a perturbation expansion about this 
state. This is the normalisation that we have chosen for equation (1) above, defining 
the Hamiltonian. In this normalisation, the partition function is often denoted A. On 
the simple cubic lattice, the Ising model expansion begins 

A = 1 + u^// + 3wV^ + • • • (7) 

As noted above, we express the field variable as /j, = 1 — x and, in order to reduce the 
memory required by the computer program, truncate the field dependence at so 
that 

A = Ao (u) + xAi (u) + x^K2 {u) + ... (8) 

We express Aq as 

oo 

Ao = ^A„«" (9) 

and define an order parameter 



oo 



M = l + ^^ = Vm,«" (10) 

^ ~ ^ ^0 n=0 



and susceptibility 



A2 Ki [KiY 

We have determined the low-temperature expansion using cross-sections up to 4 x 5, 
giving Ao, M and x correct to u^^ . The coefficients A„, mn, Cn are listed in Table 2. 

The calculations build up the finite lattice one site at a time, running through 
all q states of the site and applying a weight of 1 for state '0' and ji (expressed as 
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1 — x) for all other states. The bonds linking each new site to the partly constructed 
lattice are given weight 1 if the sites at each end are in the same state and weight 
z otherwise. For low-temperature expansions, the set of bonds considered includes 
bonds connecting the finite lattice to an outer boundary of sites taken as being in 
state '0'. 

The high-temperature expansions have been obtained using cross-sections up to 
4x5, giving Z to order v'^'^. While table 1 indicates that the vector size increases by a 
factor of 16 (or more generally q^) on going from 4 x 4 to 4 x 5, for high temperature 
expansions, the memory requirements are reduced by a factor of 2 (or more generally 
a factor of q) by making full use of the symmetry and by a further factor of 3 since the 
high-temperature field dependence is truncated at zeroth order compared to the 
truncation of the low-temperature series. In principle it is possible to further reduce 
the memory requirements for the high-temperature Ising free energy calculation by 
making use of the ferromagnetic/antiferromagnetic symmetry of the zero-field Ising 
model to reduce the vector sizes by a further factor of 2. Since we have been using 
programs designed for general q, we have not exploited this possibility. 

For the high-temperature expansions, the field weighting is not included and the 
set of bonds does not include any bonds extending beyond the finite latttice. However 
the main difference from the low-temperature expansion arises from the relation 
z = (1 — v)/{l + {q — l)v). For a finite lattice a, rather than expand Aq, we expand: 

$a=g-^(")(l + (g-l)^;)^(")A« (12) 

whence 

= g.(a)-b(a)(l + _ 1)2)K«) (13) 

where b{a) and s{a) are the number of bonds and sites respectively in lattice a. The 
expansion of $a is obtained by giving bonds a weight of 1 -|- (g — l)v for pairs of sites 
in the same state and weight 1 — v otherwise. An additional weight of q~^ is applied 
at each site. 

For the infinite lattice limit, we have 

q^-''/^{l + {q-l)zy/^^v) = q{l + {q-l)v)-''/^^v) = A^ "Q ^y(a) ^^4) 

aeB{s) 

where f is the lattice coordination number so that v/2\s the number of bonds per site, 
which is 3 on the simple cubic lattice. By using the facts that X]aGB(s) b{oL)y{oL) = ^/'^ 
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and '^aeB{s) s{a)V{a) = 1 we obtain the relation 

Hv)^ n ^o.ivv^''^ (15) 

ae-B(s) 

It is this last expression that we use in our high-temperature calculations. 

For the Ising model, it is usual to shift the zero of energy so that parallel and 
antiparallel pairs have energies ±J with J = ^^E. This leads to the more familiar 
form of the Ising model expansion: 

Z = expi+J/kT)"/^ A = 2[cosh( J/A;T)]^/2 ^^^^ ^^g) 
The high-temperature expansion for the free energy is written as 

^v) = J2(^nv'' = l + 3v^ + ... (17) 

n 

(or more generally l + 3{q — l)v'^ . . .). The coefficients a„, for n < 22 are listed in Table 
2. 

The internal energy, U, is derived from the free energy, F, or the partition function 
A by 

d d 

For low temperatures we use 

d _ dz d _ d 
dP^dpd^^~ ^di 

to give 

U = AEz^/A 
dz 

For high temperatures we use z = {l — v)/{l + {q — l)v) and ^ = —q/ (l + (? — 1)^)^ 
to give 

d ^^(l- '')(l + (^i-l)'') d 



dp q dv 

whence 

C/ = AE^i^(l - ^) - A^(l-^)(l + (g-l)^) d in^(^) 
2 q ^ ' q dv ^ ' 

The high-temperature and low-temperature expansions were computed using two 

very similar Fortran programs written for general (integer) q. Our programs work in 

terms of the general variables 2; or v, even in those special cases where only even powers 
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occur. Most of the other technical aspects are similar to our previous calculations 
using the finite lattice method. The calculations were performed using arithmetic 
modulo various primes, slightly below 2^^. Calculations were performed using 32- 
bit integers and the large vectors of residues stored as 16-bit integers. The factor 
[q — required in the general form of M is calculated as {q — l)f~^ modulo p 

by virtue of Fermat's theorem. Similarly the factors q~^ required in evaluating the 
$a were expressed as modulo p. Since all of the $q, and Aq, are of the form 

1 + kv^ -\- . . . and 1 -|- (g — l)^z^ + . . . respectively, the expansions for negative powers 
of the $a and Aq, can readily be calculated. The programs were run on an IBM 
3090/400 J with ^ GB of main memory and 2 GB of extended storage. 

4. Analysis 

Before presenting the analysis of our new series, we note a number of relevant prior 
studies. Gaunt and Sykes (1973) extended the diamond and f.c.c. low-temperature 
series, and obtained 0.307 < /? < 0.317 and a value for 7 around 1.27 - 1.30. They 
noted that the sequences of exponent estimates were converging rather slowly. A range 
of previous work was consistent with these estimates. At high temperatures, Sykes 
et ai (1972) and Camp et ai (1976) found estimates for a consistent with 0.125, in 
accordance with the belief at the time that a was a simple fraction, probably |. 

The field theory estimates of Le Guillou and Zinn- Justin (1980), in which scaling 
and high/low-temperature exponent symmetry are implicit, are 7 = 1.241(2), a = 
0.110(5), f3 = 0.325(2). Recently Oitmaa et al (1991) obtained extended low- 
temperature series for the (2-|-l) dimensional Ising model, equivalent to a Z2 gauge 
model in (2-1-1) dimensions, which is in the three-dimensional Ising universality class. 
They also derived and analysed series on the triangular and honeycomb lattices. 
Unbiased exponent estimates were given as 7 = 1.28(2), /3 = 0.311(4) and a = 0.11(4), 
while biased estimates were given as 7 = 1.25(2), /3 = 0.318(4) and a = 0.096(6). 

Nickel (1991) recently re-examined the field theory exponents, and showed that, 
by permitting a second confiuent exponent, critical exponent estimates were obtained 
that agreed rather well with those obtained from high-temperature series expansions. 
Nickel's preferred values are 7 = 1.238, u = 0.630, rj = 0.0355. Hence a = 0.110 and 
P = 0.326. The critical exponent Ai = (jjiu = 0.53. 
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Recently Mojumder (1991) has developed a theory based on partial non- 
renormalization of superconformal dimensions of matter fields on a (2,0) 
supersymmetric string world sheet. A consequence of this theory is the predicition of 
critical exponents which are supposed to be exact. These are cu = |, /3 = |, from 
which follows 7=1 and = §• 

Before analysing the newly obtained series, we remark that the analysis of the low- 
temperature series is made more difficult than the analysis of their high-temperature 
counterparts by the presence of non-physical singularities closer to the origin than 
the physical singularity. While differential approximants effectively afford analytic 
continuation beyond the nearest singularity, the series coefficients are nevertheless 
dominated by the non-physical singularity. Another approach to this difficulty is 
to apply a transformation, to move the physical singularity closer to the origin 
than the non-physical singularity. However, such 'singularity- moving' transformations 
introduce long-period oscillations (Hunter, 1987; Guttmann 1989) thus rendering 
suspect extrapolations based on such transformed series as in the analysis by Bhanot 
et al (1992). 

Another problem with the analysis of the three-dimensional Ising series is the 
presence of confluent singularities (which are extremely weak or non-existent in the 
two-dimensional model). Both field theory and high-temperature series analysis 
suggests a value of the confluent exponent not very different from 0.5. Roskies (1981) 
introduced an effective transformation for analysing such series, in which one replaces 
the original expansion variable ^ by a new variable y, defined by 1 — y = (1 — z/ Zc)^/'^ . 
If the original series had square root correction terms, the transformed series has 
analytic correction terms. If there were no square root correction terms, then nothing 
is lost by applying this transformation. One problem is that the critical temperature 
must be accurately known, but for the three-dimensional Ising model this is the 
case. The consensus of extensive high-temperature series and Monte Carlo work is 
tanh( J//cTc) = 0.218093, or Uc = exp(-4J//cTc) = 0.4120494, with errors being a few 
parts in the sixth decimal place (Guttmann 1987b). 

Our initial analysis of the three low-temperature series (specific heat, 
magnetisation and susceptibility) by Pade approximants and differential approximants 
was not particularly illuminating. Exponent estimates were a' = —0.1 to 0.2, = 0.30 
to 0.32 and 7' = 1.2 to 1.3. Biased estimates (with the critical point specified) were 
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a' ^ 0.20, P ^ 0.32 and 7' ^ 1.25. 

Analysing the Roskies-transformed series by evaluating Dlog Fade approximants at 
y = 1 gave the results shown in Tables 3 to 5. The approximants are fairly consistent, a 
few wildly deviant ones being due to defective approximants. For all three exponents 
we have taken the mean of the last few values (from 14 to 16 approximants), and 
ignored the wildly different outliers. The mean of these approximants, with error given 
as one standard deviation, is a' = 0.124 ±0.006, /3 = 0.329 ±0.009, 7' = 1.251 ±0.028. 

The effect of allowing for the correction-to-scaling shows that its neglect in the 
analysis by Bhanot et al (1992) is an even more serious problem for their analysis 
than the use of the Euler transformation. 

It can be seen that, with the exception of the estimate of a', the exponent estimates 
above are consistent with the field-theory estimates. The value for a' just fails to 
overlap the range of the field theory estimate, but we believe that the high-temperature 
series result (see below) must be taken into account also. 

We have also tested the sensitivity of our estimates to changes in both the 
critical point estimate, and the exponent estimate used in the Roskies transformation. 
Changing the estimate of the critical point by 2 parts in the fifth significant digit 
produced a change in the exponent estimate of one-tenth of the error estimates quoted 
above. That change in the critical point is far greater than the uncertainty in the 
critical point, so that we can safely assert that the results quoted include errors due 
to the uncertainty in the critical point. Changing the exponent estimate to 0.53 from 
the value 0.5 used in the Roskies transformation produced an error that was less than 
one-fifth of the errors quoted above, so that the uncertainty in the correction-to-scaling 
exponent is also negligible in this analysis. 

We have also looked carefully at the non-physical singularity that is closer to the 
origin than the physical singularity. We find its position to be at tt = — 0.2853±0.0003, 
and the exponents of the susceptibilty, magnetisation and specific heat to be 1.03, 
—0.01 and 1.01 respectively. It is difficult not to suggest that these are 1, and 1 
respectively, with presumably logarithmic corrections. These values also satisfy the 
same scaling relation a' + 2/3 + 7' = 2 that the exponents at the critical point satisfy, 
though the thermodynamic argument that leads to this result at the critical point is 
not obviously applicable in this non-physical region. 

For the high-temperature series, the critical temperature is the radius of 
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convergence of the series. It follows that the ratio method and its variants can be 
used. Unbiased differential approximants, using the method of analysis described in 
Guttmann (1987a) give v"^ = 0.0475185, a = 0.126. Utilising the linear dependence 
between the estimates of Vc and a gives a biased estimate of 0.11 for cu at the correct 
vl = 0.047565. 

With a 'correction-to-scaling' exponent close to 0.5, the ratio of succesive terms 
in the specific heat expansion is expected to behave ctS '^fx — /-^ [1 + (a - l)/n + 
cjn^l'^ _l_ 0{\jyi^y^ so estimators of a are given by the sequence (r„/// — l)n + 1 = 
CKi -I- c/n}/'^ + 0{l/n) where ji = In Fig. 1 we show the plot of this sequence 

against If-n}^"^. The series we have used is the expansion of ^ in powers of v^, just 
as was used by Sykes (1972), but with three further terms. Linear regression to the 
data gives a = 0.113 — 0.0637/n-'^/^-|-0.389/n. As we have shown in previous work, the 
method of differential approximants can accurately predict the most significant digits 
of the next term in the series (see e.g. Guttmann and Enting 1988a). In this way 
we have estimated the coefficient of f in the zero-field high-temperature partition 
function to be 2.73376 * 10^^, where only the last quoted digit is uncertain. Using 
the additional term in the analysis just described, we obtain an estimate for a of 
0.111, with the correction terms given by — 0.0511/n-^/^ -|- 0.364/n. This estimate of 
a is impervious to the uncertainty in the estimated coeflBcient. We conclude from 
this analysis that a = 0.110 ± 0.005. The comparatively small amplitude of the term 
proportional to helps us to understand why the series is so well behaved. 

Recently, Liu and Fisher (1990) studied the correction-to-scaling amplitudes of 
the three-dimensional Ising model and argued that the amplitudes of the correction- 
to-scaling term should be negative. Writing C ~ A|t|~"[l + ae|t|^ + ait + . . .] where 
t = (T — Tc)/Tc and ^ = | our result transforms to ae ~ —0.040, which has the 
predicted sign. 

Biased first-order differential approximants give a = 0.104 ± 0.018, where the 
analysis is precisely that described in Guttmann (1987a). The error bars correspond 
to two standard deviations. A Roskies-transformed Fade analysis was also performed 
for the high-temperature series, paralleling the low-temperature investigation. The 
short series meant that very few approximants were obtained, and so the results are 
not particularly helpful (exponent estimates in the range 0.08 to 0.3 were obtained). 
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5. Concluding remarks 

The calculations presented here have shown that on the simple cubic lattice 
the finite lattice method of series expansion compares favourably with conventional 
expansion methods. However there is not the dramatic difference that occurs for 
many two-dimensional problems. Nevertheless, significant extensions of both high- 
and low-temperature series have been obtained. 

Our exponent estimates are all consistent with field theory exponent estimates. 
The two best behaved series give f3 = 0.329 ± 0.009 and a = 0.110 ± 0.005. The 
estimates of a' and 7', while consistent with field theory estimates, had rather wide 
error bars, and so the agreement with field theory was less convinicing. Our estimates 
of a and /? are not consistent with the conjectured exact values of Mojumder. 

After acceptance of this paper, we received a preprint from Vohwinkel (1992), as 
well as some additional series. Vohwinkel has used a version of the shadow-lattice 
method of Sykes, and has obtained six further low temperature series coefficients. 
Analysis of these longer series by the methods used in this work do not significantly 
change our exponent estimates. 
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Tables and table captions 



Cross-section 


s 


Order of Potts (Isin{ 
Low-T 


5) series 
High-T 


R{d X w,2) 


1 X 1 


4 


13 (6) 


3(2) 


2 


1x2 


5 


17 (8) 


5(4) 


4 


2x2 


7 


25 (12) 


9(8) 


16 


2x3 


8 


29 (14) 


11 (10) 


64 


3x3 


10 


37 (18) 


15 (14) 


512 


3x4 


11 


41 (20) 


17 (16) 


4096 


4x4 


13 


49 (24) 


21 (20) 


65536 


4x5 


14 


53 (26) 


23 (22) 


1048576 


5x5 


16 


61 (30) 


27 (26) 


33554432 



Table 1. Combinatorial factors determining the computational complexity of finite 
lattice series expansions on the simple cubic lattice. Column 1 gives the various 
possible cutoff points as specified by the largest cross-section that need be considered 
after making optimal use of cubic symmetry. The second column is s, the largest 
value of the sum of length plus width plus depth for any of the cuboids used in the 
expansion. The low-temperature series for the Potts model can be obtained to order 
^4s-3 ^Yie high-temperature series to v'^^~^. Ising model series can be obtained 
to tt^*~^ and The final column, i?(.,2), gives the number of vector elements 

required for the low-temperature Ising expansion. 
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n 








On 





1 


1 





1 


3 


1 


-2 


1 





4 











3 


5 


3 


-12 


12 





6 


-3 


14 


-14 


22 


7 


15 


-90 


135 





8 


-30 


192 


-276 


192 


9 


101 


-792 


1520 





10 


-261 


2148 


-4056 


2046 


11 


807 


-7716 


17778 





12 


-2308 


23262 


-54392 


24853 


13 


7065 


-79512 


213522 





14 


-21171 


252054 


-700362 


329334 


15 


65337 


-846628 


2601674 





16 


-200934 


2753520 


-8836812 


4649601 


17 


627249 


-9205800 


31925046 





18 


-1962034 


30371124 


-110323056 


68884356 


19 


6192066 


-101585544 


393008712 





20 


-19610346 


338095596 


-1369533048 


1059830112 


21 


62482527 


-1133491188 


4844047090 





22 


-199807110 


3794908752 


-16947396000 


16809862992 


23 


641837193 


-12758932158 


59723296431 





24 


-2068695927 


42903505030 


-209328634116 


27337xxxxxxx 


25 


6691611633 


-144655483440 


736260986208 





26 


-21710041944 


488092130664 


-2582605180212 





Table 2. Coefficients in low-temperature expansions for Aq, M and x tiigfi- 
temperature expansion for $, defined by equations (9), (10), (11) and (17). Tlie 
incomplete coefficient, a24, was not derived using the finite lattice method, but 
represents an extrapolation obtained using differential approximants. 
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N 


[N - 2/N] 


[N -1/N] 


[N/N] 


[N + 1/N] 


[N + 2/N] 


5 


0.4192 


0.3590 


0.2290 


0.5032 


0.3240 


6 


0.3261 


0.3346 


0.3431 


0.3437 


0.3325 


7 


0.4559 


0.3437 


0.3431 


0.3332 


0.3326 


8 


0.3027 


0.3388 


0.3258 


0.3342 


0.3322 


9 


0.3358 


0.3419 


0.3731 


0.2900 


0.3279 


10 


0.3275 


0.3186 


0.3580 


0.3215 


0.3230 


11 


0.3270 


0.3318 


0.3233 


0.3210 


0.3391 


12 


0.3290 


0.3274 


0.2624 


0.2113 


0.3286 


13 


0.3286 


0.2154 









Table 3. Estimates of /3 using Dlog Pade approximants to Roskies-transformed series. 



N 


[N - 2/N] 


[AT - 1/iV] 


[N/N] 


[N + 1/N] 


[N + 2/N] 


5 


1.1962 


1.2342 


1.2997 


1.3836 


1.1646 


6 


1.4124 


1.4975 


1.3363 


1.2835 


1.1892 


7 


1.4200 


1.3068 


1.3431 


2.5047 


1.1088 


8 


1.7242 


1.2416 


1.2809 


1.2551 


1.2049 


9 


1.2981 


1.2655 


1.2940 


1.2087 


1.2047 


10 


2.7628 


1.2466 


1.2493 


1.2767 


1.2511 


11 


1.2494 


1.2461 


1.2605 






12 


1.2643 











Table 4. Estimates of 7' using Dlog Pade approximants to Roskies-transformed series. 
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N 


[N - 2/N] 


[N -1/N] 


5 


0.1850 


0.1197 


6 


0.1078 


0.1173 


7 


0.1224 


0.1131 


« 


n 1260 


n 1207 


9 


0.1234 


0.1686 


10 


0.1221 


0.1214 


11 


0.1221 


0.1182 


12 


0.1375 





Table 5. Estimates of a' using Dlog 
series. 



[N/N] 


[N + 1/N] 


[N + 2/N] 


0.1102 


0.1410 


0.0829 


0.1236 


0.1319 


0.1265 


0.1284 


1329 


1 109 


0.1239 


0.1261 


0.1364 


0.1226 


0.1208 


0.1658 


0.1227 


-0.686 


0.2232 


0.1743 







e approximants to Roskies-transformed 
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Figure captions 



Figure 1. Ratio method estimates of a, plotted against v}^'^. The series used was 
expanded in powers of v^. The slope gives the amplitude of the 'correction-to-scaling' 
term and the intercept gives the estimated value of a after taking the correction-to- 
scaling into account. 



24 



